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ABSTRACT 

Using spectroscopic observations and photometric light curves of 238 nearby M dwarfs from the 
MEarth exoplanet transit survey, we examine the relationships between magnetic activity (quantified 
by Ha emission), rotation period, and stellar age. Previous attempts to investigate the relationship 
between magnetic activity and rotation in these stars were hampered by the limited number of M 
dwarfs with measured rotation periods (and the fact that v sin i measurements probe only rapid 
rotation). However, the photometric data from MEarth allows us to probe a wide range of rotation 
periods for hundreds of M dwarf stars (from shorter than than one to longer than 100 days). Over 
all M spectral types that we probe, we find that the presence of magnetic activity is tied to rotation, 
including for late-type, fully convective M dwarfs. We also find evidence that the fraction of late- 
type M dwarfs that are active may be higher at longer rotation periods compared to their early-type 
counterparts, with several active, late-type, slowly rotating stars present in our sample. Additionally, 
we find that all M dwarfs with rotation periods shorter than 26 days (early-type; M1-M4) and 86 
days (late-type; M5-M8) are magnetically active. This potential mismatch suggests that the physical 
mechanisms that connect stellar rotation to chromospheric heating may be different in fully convective 
stars. A kinematic analysis suggests that the magnetically active, rapidly rotating stars are consistent 
with a kinematically young population, while slow-rotators are less active or inactive and appear to 
belong to an older, dynamically heated stellar population. 

Subject headings: stars: low-mass — stars: activity — stars: late-type — stars: rotation — stars: 
kinematics — stars: chromospheres 


1. INTRODUCTION 


Many late-type M dwarfs (> M4) have strong mag¬ 
netic fields that can exceed 1 kG and heat their stellar 
chromospheres and coronae, creating “magnetic activ- 
ity” that is observe d from the radio t o the X-ray (e.g., 
Hawley e t al.l 119 9^ IWest et alJ 120 04 IRei ners &: Basril 


I 


20081 Berger et al.| I2008[ IReiners fc Basril 1200 

Williams et al.l 120141) . Although magnetic activity 
has been observed in M dwarfs for decades, the exact 
mechanism that gives rise to the chromospheric and 
coronal heating is still not well-understood. The produc¬ 
tion of magnetic fields and the subsequent activity may 
play a vital role in the habitability of attending planets, 
which may be numerous given the u biquity of M dwarfs 
as pl a net hosts in the Galax y fe.g.. iCh arbonneau et ah 
2009 : iMuirhead et al.l 120121 : iDressing &: Charbonneau 
2013 . 1201511 . In addition, given their fully convective 
interiors, late-type M dwarfs serve as important lab¬ 
oratories for studying magnetic dynamo generation 
in stellar (and potentially planetary) environments 
and are vital for understanding the role that various 
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stellar properties play in the magnetic field and activity 
generation in low-mass stars. 

In solar-type stars, magnetic field generation and 
subsequent heating is closely tied to stellar rota¬ 
tion; the faster a star rotates, the stronger its mag¬ 
netic heating and surface activity. Angular mo¬ 
mentum loss from magnetized winds slows rotation 
in solar-type stars, and as a result, magnetic ac¬ 
tivity decreases with age. These effects have been 
studied for decades, producing ample evidence for a 
strong connection between age, stellar ro tation, and 
magnetic activity in solar-type stars (e. g .. iSkumanichl 
197^ lBarrvlll98^ iSoderblom et all 119911: iBarnesI 20031 : 


Pizzolato et al.l 120031 : iMamaiek fc Hillenbrandl 120081) . 

All indications suggest that this connection between age, 
rotation and activity extends to early-type M dwarfs 
(< M4), w here rotation and activity are strongly corre¬ 
lated (e.g..lPizzola.to et al.ll2003l : iMohantv fc Basrill2003l : 
iKiraea fc Steoienl l2007f) . The finite active lifetimes of 
early-type M dwarfs observed in nearby clusters suggest 
that age continues to play an important role in the ro- 
tation and magnetic activity evolution o f low-mass stars 
(|Stauffer et al.l r994 iHawlev et al.lll999ll . 

At a spectral type of about M4, stars become 
fully convective (0.3 5 Mq; iGhabrier fc Bara^ 119971 : 
iReid fc HawlevI [20051 1. a property that may affect how 
magnetic field (and the resulting heating) is generated. 
Despite this change, magnetic activity persists in late- 
type M dwarfs; the fraction of active M dwarfs peaks 
around a spectral type of M7 before decreasing into 
the b r own dwarf regime (|Hawlev et al.l[T9^ iGizis et ‘ahl 
120001 : IWest et HI 1200411 . Several studies have demon- 
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strated that the large fraction of late-type M dwarfs ob¬ 
served to be active is likely a r esult of their long activity 
lifetimes (|Silvestri et alJ[2005l : iWest et al.ll2008ii . 

Using 37,845 M dwarfs from the Sloan Digital Sk y Sur- 
vey (SDSS; I York et al .1 [2nn(ill . IWest et aiJ jlOOl, iMol) 
found that stars farther from the Galactic plane were less 
likely to be magnetically active (as traced by Ha) than 
those near the Plane. They interpreted the change in ac¬ 
tivity fraction as an effect of age: as stars pass through 
the Galactic plane, they are dynamically heated by grav¬ 
itational interactions and obtain orbits that take them 
farther from the Plane. Therefore, stars close to the 
Galactic plane are statistically younger, while sta rs far¬ 
ther away are statistically older. IWest et al.l 1)200811 quan¬ 
tified this “Galactic stratigraphy” using a ID dynamical 
model and derived the Ha activity lifetimes for M0-M7 
dwarfs, finding that early-type M dwarfs (with both ra¬ 
diative and convection zones) have active lifetimes of 1-2 
Gyr, while late-type M dwarfs (with fully convective in¬ 
teriors) have active lifetimes that exceed 7 Gyr. The level 
of activity (as quantified by the ratio of the luminosity in 
Ha to the bolometric luminosity — Tna/Tboi) also ap¬ 
pears to decrease as a function of stratigraphic age for all 
M dwarfs (early and late-type), confirming that an age- 
activity relation p ersists into the fully-convective regime 
(IWest et al.ll2n08ll . 

Tying activity (and age) to rotation for late-type M 
dwarfs has been more challenging. Recent simulations 
of dynamos find that rotation may play a significant 
role in the magnetic held generation of fully convec¬ 
tive stars (|Dobler et al.l 120061: iBrownind 1200^. Froin 
a simple analytical model, iReiners fc Mohantvl (|2012ll 
suggest that the age-rotation relation extends to late- 
type M dwarfs and that the angular momentum evolu¬ 
tion of all stars is more a function of stellar radius than 
it is interior structure. Indeed, a few empirical stud¬ 
ies have uncovered evidence that activit y and rotation 
migh t be linked in late- t ype M dwarfs (iDelfos se et al.l 
19981: iMohantv fc Basril 2003t IReiners fc Basril 120081 : 


Browning et al 


V &: jjasril UUUdt IKeiners jjasril izuust 
.1 120101 : IReiners fc Mohantvl l2012fl . How¬ 


ever, the majority of the observations have relied on 
high-resolution spectroscopic data, which measure the 
rotational velocity modified by the inclination of the star 
(y sin i). While v sin i measurements provide important 
clues to the underlying stellar rotation, they have a cou¬ 
ple of limitations: 1) the best spectrographs can produce 
V sin i values only down to about 1 km s“^, which for a 
0.2 Rq star corresponds to a rotation period of 10 days, 
and are therefore completely insensitive to slowly rotat¬ 
ing stars; and 2) the inclination dependence adds signif¬ 
icant scatter to the derived rotation velocities. Previous 
results demonstrate these limitations by showing that al¬ 
most all late-type M dwarfs with a detected v sin i have 
the same “saturated” level of magnetic activity; slowly 
rotating M dwarfs with potentially less magnetic activity 
are not detectable. 

An alternative method for studying stellar rotation 
is to use photometricall y derived rotation periods (e.g., 
iKiraga fc: Steoienl l2007tl . Periodic signals result from 
brightness variations caused by long-lived spots on the 
stellar surface rotating in and out of view. While v sin i 
observations require a single spectroscopic observation, 
observing periodic, photometric variability requires nu¬ 


merous, high-cadence observations that have been histor¬ 
ically prohibitive for large samples of late-type M dwarfs. 
However, recent programs to search for transiting planets 
around late-type M dwarfs have produced large catalogs 
of time-domain photometry from which can be gleaned 
several i mportant stellar properties, includ ing rotation 
periods (|Irwin et al.ll2009a IlW et al.ll2012ll . 

One of these transit pro g rams, t he MEarth P roieclLl 
(iNutzrn an fc Gharbon nead 120081 : lirwin et al.l I2009H 
I2011bl: IBerta et al.l I2012r is employing two arrays of 
robotic telescopes to photometrically monitor 4,000 
nearby, mid-to-late M dwarfs. MEarth data and as¬ 
sociated observations have provided new measurements 
of the fundamental properties of nearby mid-to-late M 
dwarfs, including their distances (lDittnianiL_eL_al 1201^ 
and their near-inf rared spectra iNewton et al.l 201^ . 
Early in the survey, lirwin et al.l (|2011afl determined pho¬ 
tometric rotation periods (ranging from 0.28 to 154 days) 
for 41 of the MEarth M dwarfs with known parallaxes, 
and identified a relationship between these rotation pe¬ 
riods and kinematic age. As the mid-to-late M dwarfs 
observed by MEarth are nearby and therefore relatively 
bright, it is possible to obtain a larger sample of optical 
spectroscopy to probe magnetic activity, thus enabling a 
robust investigation of how the rotation-activity relation 
for fully convective M dwarfs extends to slowly rotating 
stars. 

Determining ages for stars can be particularly cha l- 
lenging, especially for low-mass stars (|Soderblomll2010ll . 
In the past, age-dependent relations for stars have been 
calibrated using stellar clusters. However, the faintness 
of M dwarfs and the large distances to stellar clusters 
with ages > 1 Gyr preclude any detailed observations 
of low-mass stars in the cluster environment. One al¬ 
ternative method for probing stellar age is to use the 
3D stellar kinematics of a population. Historically, the 
dynamics of stars have been used to estimate ages, as 
dynamical interactions cause the orbits of older stars to 
become more elliptical an d inclined to the Galactic plane 
(iWiEon fc Wooll^ 119701 : lEggen fc IbenI 119891 : iLeggettI 
Il992f) . Specihcally, stars that exhibit larger velocity dis¬ 
persions have likely experienced more dynamical encoun¬ 
ters with molecular clouds and/or other stars and are 
therefore older. This is particularly true in the vertical 
or W direction of stellar motion since most stars begin 
their lives in the Galact ic disk and slowl y diffuse away 
via dynamical heating (|West et al.l 1200611 . Older stars 
also have slower azim uthal or V mean vel ocities due to 
asymmetric drift ('e.g.. lEggen fc Ibe'i]ll989H . While using 
stellar kinematics to determine ages of individual stars is 
highly problematic (since the dynamical processes that 
affect individual stellar orbits are stochastic), 3D mo¬ 
tions can reveal important age information for bulk pop¬ 
ulations. Dynamical analyses do require large samples, 
which until recently have not been available for late-type 
M dwarfs. 

In this paper, we use an expanded sample of MEarth 
rotation periods for nearby M dwarfs, new optical spec¬ 
troscopy of these targets, and their full space motions 
to investigate the rotation-activity and rotation-activity- 
age relation for late-type M dwarfs. Section 2 describes 
the determination of rotation periods from MEarth data 

^ http://cfa.harvard.edu/MEarth/ 
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and the determination of magnetic activity from optical 
spectroscopic observations. We describe the rotation, ac¬ 
tivity and kinematic analysis in Section 3. In Section 4 
we demonstrate that there is a close tie between rotation 
and activity and that the slowest rotators come from a 
dynamically older population. Lastly, Section 5 summa¬ 
rizes our conclusions and gives a brief discussion. 

2. DATA 

2.1. Rotation Periods from MEarth Photometry 

To measure M dwarf rotation periods, we use photo¬ 
metric observations fr om the MEarth survey for tran¬ 
siting exoplanets (see iNutzman k, Charbonneaul 120081 : 
iBerta et al.ll201^ . MEarth comprises two telescope ar¬ 
rays, each consisting of eight robotic telescopes. All of 
the data in this paper come from the MEarth-North ar¬ 
ray at the Fred Lawrence Whipple Observatory (FLWO) 
at Mt. Hopkins, AZ, which has been gathering data 
since 2008. The MEarth team recently commissioned 
a duplicate array at Cerro Tololo Interamerican Obser¬ 
vatory (CTIO), called MEarth-South. Because the new 
array started regular observations only in January 2014, 
data from MEarth-South are not included in the follow¬ 
ing analyses. The MEarth M dwarf light curves used in 
this paper (along with the full MEarth sample) are pub¬ 
licly available for download from the MEarth websitcQ. 

The sample of M dwarfs observed by MEarth was 
designed to maximize sensitivity to the detection of 
small transiting planets (2 — dRm), as detailed in 
INutzman fc Charbonneaul (|2008ll . The MEarth sample 
is heavily biased toward mid-to-late M dwarfs within 
33pc, with most stars having spectral types of M4-M6 
but also including several stars of earlier and later spec¬ 
tral types. This paper makes use of the MEarth sample 
of stars and inherits the selection criteria of the planet- 
search survey. While the sample of stars observed was 
not designed explicitly to probe rotation across the en¬ 
tire M dwarf spectral class, it spans a sufficient range of 
spectral types to enable interesting comparisons between 
earlier (which in this paper refers primarily to M3-M4) 
and later (primarily M5-M6) spectral types. 

The MEarth telescopes are 40 cm in diameter and 
equipped with back-illuminated CCDs. The telescopes 
observe automatically whenever conditions allow, with 
an adaptive scheduling algorithm selecting targets from 
a prioritized queue. The rotation periods in this paper 
are derived from “planet cadence” observations, in which 
a star is observed once every 10 — 30 minutes for at least 
one observational season, resulting in light curves con¬ 
taining from hundreds up to a few thousand data points. 
For some stars, up to three observational seasons con¬ 
tributed to the rotation analysis. We took care to min¬ 
imize changes to the telescopes to preserve long term 
stability, b ut those changes th at were necessary are sum¬ 
marized in IBerta et al.l (|2012f) and documented in detail 
in the notes associated with the MEarth data releases. 

The MEarth-North bandpass is fixed and nearly iden¬ 
tical across the eight telescopes. Designed to maximize 
photon flux from red M dwarf targets, the bandpass is 
wide and has a shape set at short wavelengths by a long- 
pass hlter that cuts on at 715 nm and at long wave¬ 
lengths by the quantum efficiency of the CCD, which 


decreases significantly by about 1000 nm. The time- 
variable optical depth through strong telluric water vapor 
features that fall within the MEarth bandpass can cause 
systematic photometric trends in MEarth M dwarf light 
curves, but the ef fect can be mitiga ted through an em¬ 
pirical correction (jBerta et al.ll20l3l . Photometric mod¬ 
ulations due to rotating starspots are relatively muted 
in the MEarth bandpass compared to bluer wavelengths, 
but they can still have readily dete ctable peak-to-pea k 
amplitudes of up to several percent (|Irwin et al.ll20il^ . 

Rotation periods were extracted using a weighted least- 
squares periodogram fitting algorithm that corrects for 
and marginalizes over systematic photometric trends, in¬ 
cluding those caused by varying telluric w ater absorp¬ 
tion. T his analysi s is described in de tail in lirwin et al.l 
(|2011ajl . Whereas H rwin et rdl (j2011a[l included rotation 
periods only for stars with published literature paral¬ 
laxes, here we extend the sample to all M dwarfs with 
rotational modulations detectable in the MEarth pho¬ 
tometry that were available up to 2011, resulting in 164 
targets. For this paper’s analysis, we applied a two-tier 
classification system (“1” and “2”) for all of the measured 
rotation periods, with “1” being the most robust detec¬ 
tions. Some of the reasons that rotation periods were 
given a classification of “2” were if the stars had long pe¬ 
riods and less than two complete cycles were observed or 
if the stars had larger uncertainties in the phase folding 
procedure. We were unable to measure rotation periods 
for all of the MEarth stars. For stars with detected ro¬ 
tation period, the amplitude of the photomet ric modula¬ 
tions w ere typically 0.5 — 2% peak-to-peak. lirwin et al.l 
(l2011al) made initial attempts to characterize the reli¬ 
ability and contamination of rotation period detection 
with MEarth and found the selection to be fairly clean 
for stars with >1% variability, but messy (and strongly 
period-dependent) for stars showing variability of 0.5% 
or less. Figure [T] shows the distribution of measured rota¬ 
tion periods for all of the M dwarfs with spectroscopic ob¬ 
servations. The periods used in this paper can be found 
in Table [T] A complete catalog of rotation periods from 
the MEarth survey will be included in a future publica¬ 
tion (Newton et al. in prep). 



log[Period] (Days) 

Fig. 1.— Distribution of rotation periods for all MEarth M 
dwarfs (solid) and late-type (> M4) M dwarfs (dashed) with spec¬ 
troscopic observations. 
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2.2. Magnetic Activity from FAST Spectroscopy 

To probe the magnetic activity of this sample, we ob¬ 
tained optical spectra for 238 MEarth M dwarfs. The 
sample of stars observed spectroscopically included those 
with measured rotation periods, as well as additional 
M dwarfs for which no rotational modulations were de¬ 
tected. Observed using the 600 lines/mm grating on the 
FAST spectrograph on the 1.5m Tillinghast Telescope at 
FLWO, the spectra have a resolution of i? = 3000. The 
spectra cover 5550-7550 A and include several features 
that are important for M dwarf analyses, including (but 
not limited to) the CaH and TiO molecular bands and 
the Ha atomic line. 

The FAST observations were acquired over 30 nights 
from December 2010 to July 2012, with HeNeAr lamp 
exposures taken at every telescope position and spec- 
trophotometric standards observed every night. Expo¬ 
sures lasted typically about five minutes, yielding signal- 
to-noise ratios (SNR) of about 50 per resolution element. 
With a scripted IRAF reduction, we bias-subtracted and 
flat-fielded all spectra using calibration exposures taken 
every time the instrument grating changed, and applied 
a wavelength calibration determined for each spectrum 
from its matching HeNeAr exposure. We extracted one¬ 
dimensional spectra for each target, weighting by the pro¬ 
file in the cross-dispersion direction and using linear in¬ 
terpolation to subtract the sky. We determined a rough 
flux calibration for each spectrum using the nightly spec- 
trophotometric standards. 


100 

80 

^ 60 
<D 

E 

2 40 

20 

0 

Fig. 2.— Distribution of spectral types for MEarth stars with 
FAST spectroscopic observations, including stars with measured 
rotation periods (dashed). 
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11 stars with previ ously determined RVs from hi^h- 
resolu t ion spectra (iDelfosse et al.l 119981: iMontes et akl 
200H iNidever et al l2002t Mohantv &: Basril 120031 : 


Browning et al.l 12010 ; iShkolnik et al.l 1201011 . The RMS 
difference between the RV measurements was 5.1 km s“^, 
which we took as the uncertainty in our ability to mea¬ 
sure RVs from the FAST spectra. Each spectrum was 
corrected to air wavelengths in a zero-velocity heliocen¬ 
tric rest frame for the measurements of spectral features 
(see below). 


2.3. Distances from MEarth Astrometry and Other 
Sources 

To characterize the basic properties of the target M 
dwarfs, we used distanc e mea surements from several 
sources. iDittmann et al.l (|2014ll used MEarth imaging 
to measure geometric parallaxes for 1,507 M dwarfs in 
the MEarth sample. These include 150 out of the 164 
stars in the rotation sample and 213 out of the 238 M 
dwarfs in the spectroscopic sample. For the remaining 
stars witho ut MEarth parallax es, we adopted distances 
published in lLepine et al.l (|2005ll . which come either from 
literature trigonometric parallaxes or from photomet¬ 
ric parallax estimates. Proper motions for all of the 
star s were measured as p art of the LSPM-North cata¬ 
log (ILenine k. Sharal[200^ . 

3. ANALYSIS 

The FAST spectra w ere processed with the Hammer 
spectral-typing facility (iCovev et al.ll2007ll . The Ham¬ 
mer uses measurements of atomic and molecular features 
to estimate an initial spectral type. We then visually in¬ 
spected all of the spectra with the Hammer (v. 1_2_5) 
“eye check” function, and manually assigned spectral 
types to each star -- we identified and classified 238 M 
dwarfs, (see Figure [2]). Of the 238 M dwarfs for which 
we obtained FAST spectroscopy, 164 have measured ro¬ 
tation periods. 

We measured radial velocities (RVs) by cross- 
correlating each FAST M dwarf spectrum (between 5600- 
7100 A) with the respective iBochanski et al.l (l2007bll M 
dwar f template (|Moliantv fc Basril 120031 : iWest k. Basril 
1200911 . We determined the typical uncertainty in our 
radial velocity by comparing our calculated RVs for 


3.1. Magnetic Activity 

Using the radial velocity-corrected spectra, we deter¬ 
mined the presence and magnitude of magnetic activity 
by measuring Ha equivalent widths (EWs) as defined by 
IWest et al.1 (IMlIl . We used 6500-6550 A and 6575-6625 
A as our two continua regions, with the central Ha wave¬ 
length at 6562.8 A. In previous studies, stars with Ha 
EWs g reater than 0.75 A were considered magnetically 
active (|West et al.l 1201 111 . However, given the relatively 
small size of our FAST sample, we were able to visually 
inspect each spectrum and classified all stars with de¬ 
tectable Ha emission as “magnetically active.” All of the 
active stars in our sample have Ha EWs > than the 0.75 
A criterion from lWest et al.l (1201 If) , and would have been 
selected as “active” using an automatic EW threshold. 
The high SNR of the FAST spectra allowed us to be sen¬ 
sitive to Ha EWs far below the values measured for the 
active stars in our analysis (particularly for the M3-M6 
dwarfs), indicating that we could have detected activity 
in many of the “inactive” stars had it been present. 

M dwarfs that exhibit weak magnetic activity 
can show Ha in absorption rather than emis¬ 
sion JStauffer fc HartmannI Il986l : IWalkowicz fc HawlevI 
1200911 . As such, using Ha emission as an activity in¬ 
dicator will not include M dwarfs with the weakest levels 
of magnetic activity. Throughout this paper, we define 
active stars as those that exhibit Ha in emission and are 
therefore strongly magnetically active. We examined all 
of the inactive early-type M dwarf spectra to look for Ha 
absorption, and found no evidence for strong absorption. 
In 9 of 47 of these cases, there were hints of weak Ha 
absorption. However, a strong nearby TiO absorption 
feature precluded firm detections, given the low resolu- 
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tion of the FAST spectra. Activity fractions reported in 
this paper may be slightly higher in reality, because we 
may have excluded weakly active M dwarfs that do not 
show Ha in emission. 

We identified 160 active M dwarfs, which comprises 
67% of the M dwarfs in the spectroscopic sample. Figure 
[3] shows the fraction of magnetically active M dwarfs in 
our sample as a function of spectral type. As seen in 
previous st udies (|Hawlev et al.l Il996t iGizis et all 120001 : 
IWest et aUEiril . we see an increase in activity with in¬ 
creasing spectral type, which is likely due to the longe r 
active lifetimes of late-type M dwarfs (|West et al.ll2008[ l. 



Fig. 3.— The fraction of magnetically active M dwarfs as a func¬ 
tion of spectral type. Errors were computed from binomial statis¬ 
tics and the numbers above each symbol indicate the total number 
of stars (active -|- inactive) in each bin. Similar to previous results, 
we see a dramatic increase in the fraction of active M dwarfs with 
increasing spectral type. 


There were previously det ermined Hq EW v alues for 23 
of the stars in our sample (|Reid et al.l 1199511. W e com - 
pared the FAST EWs to those from iReid et al.l ()1995[ 1 
and found that 12 of the measurements are within 13% 
of those previously reported with the median difference 
being 14% for all stars. Six of the stars have EWs that 
are different by more than 25% and only two show vari¬ 
ations larger than 50%. Variation on this level can be 
explained b y the typical var iation of Ha emission in M 
dwarfs le.g.. lBell et al.l[2012f) . With the activity criterion 
used in this paper, none of the 23 stars switch from being 
active to being inactive (or vice versa) between the two 
epochs. 

For the active stars in the sample (those with Ha 
emission), we calculated the magnetic activity strength, 
THr>/T hni using the x factor from iWalkowicz et al.l 
(j2004f l. We multiplied the spectral type-dependent con- 
yersion factor y by the Ha EW to get the ratio THa/Tboi- 
Figure 0] shows the THa/Tboi yalues as a function of spec¬ 
tral type. The dashed line in Figure |4] indicated the 
typical level of Ha emission to which we are sensitive 
as a function of spectral type. This detection thresh¬ 
old was estimated using the typical noise level near Ha 
in the FAST spectra and assuming a detectable emis¬ 
sion line would have a peak that is at least three times 
larger than the noise. As shown in previo us studies (e.g., 
iBurgasser et al.l 120021 : IWest et all 1200411 . we see a de¬ 


crease in THa/Tboi with increasing spectral type (well 
above detection threshold for M3-M6 dwarfs; dashed 
line) as well as a large spread of THa/Tboi at each spectral 
type. The latter is a combination of the intrinsic spread 
of the active M dwarf population and age-dependent ef¬ 
fects on the activity that are washed out when stars of 
different ages are binned together. 



Spectral Type 

Fig. 4.— L^a/Lho\ a function of spectral type for active stars 
in the sample. Spectral types earlier than M3 are not shown, as 
our sample contains no such stars that are active. Median values 
of log(LHa/tjbol) are listed below. The dashed line indicates the 
typical LHa/tl/bol values to which we are sensitive as a function of 
spectral type. There is a large spread in the THa/^bol values at 
each spectral type, as well as a significant decrease in the median 
^Ho/ibol values with increasing spectral type. 


As part of our analysis, we also measured the TiO, 
CaH, and CaOH molecular bandheads as defined by 
IReid et all (|1995ll . and the Nal atomic absorption line. 
All of the values can be found in the online version of 
Table m 

3.2. Kinematics 

We combined the stellar positions, proper motions, dis¬ 
tances and radial velocities to compute the Galactic C/, 
V, W velocities for 212 of the stars in our sample^] We 
corrected for the solar motion with res pect to the Lo¬ 
cal S tandard of Rest (11, 12, 7 km s~^: ISchonrich et al.l 
1201011 and propagated our measurement errors through 
to formal uncertainties in U, H, and W using a modified 
version of the IDL routine GAL_UVW.PRD. 

There were 159 stars in our sample that had both good 
rotation measurements (flags “1” or “2” from Section 2) 
and 3D space motions. We divided these stars into four 
different bins based on the rotation velocities (<1 days, 
1-10 days, 10-100 days, and >100 days; with 58, 47, 42 
and 17 stars in each bin respectively). 

The velocity distributions of specific stellar popula¬ 
tions can be modeled by Gaussian distributions in all 
three kinematic components. Mean space velocities and 
their respective spreads and uncertainties are often calcu¬ 
lated using fits to cumulative probability dist ributions or 
simple Gaussian fits to velocity histograms (IReid et al.l 

® We used a right-handed coordinate system with U pointed 
towards the Galactic center. 
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119951: iBochanski et al.l l2007at iReiners k, Basril I2009D . 
However, these methods are computationally intensive 
and/or heavily biased by the binning strategy. The chal¬ 
lenges associated with previous methods are amplified 
when the data are sparse, which is particularly impor¬ 
tant for our small number of slow rotators. 

Alternatively, under the assumption that the underly¬ 
ing velocity distribution is Gaussian, probabilities for a 
given model mean and standard deviation can be com¬ 
puted analytically and compared to the data using a 
Bayesian approach. For a given subset of velocities, the 
joint probability distribution for the best fit mean (/r) 
and standard deviation (cr) of the population is given 
by: 


c'}; {data}) oc p({data}; {/r, <J))p{{^i, crj), (1) 

where the first term on the right represents the likeli¬ 
hood of the data and the second term on the right is 
the prior distribution of parameters, which we took to 

be uninformative^ 

The likelihood is the product of the probabilities for 
all of the data points in a subset, where the probability 
for each datum is given by the Gaussian distribution: 


p({m,s}j;{^,cr}) 


1 

V^27r(cr2+s2) 


exp 


2{a^+sf) 


( 2 ) 


where m and s are the measured velocity and uncertainty 
of a single star, p is the probability that the datum comes 
from a Gaussian distribution described by p and cr, and 
the index i runs over all stars within a bin if it contains 
at least three stars. 

We constructed a grid of p and a values and explored 
the joint probability distribution of Equation 1 for the 
four kinematic sub-samples at each of the grid points. 
From the resulting peaks of the joint probability density 
distributions, we computed p and cr values for each ro¬ 
tation/activity sub-sample. We also calculated realistic 
uncertainties from the spread of the marginalized distri¬ 
butions for each parameter. The results of this analysis 
can be found in Tableland are discussed below in more 
detail. For the fastest rotators, the data did not provide 
a good constraint on the U velocity dispersion; that value 
was excluded from Table [2l 


4. RESULTS 

4.1. Rotation vs. Activity 

From the sample of 164 M dwarfs with measured rota¬ 
tion periods, 127 were classified as magnetically active. 
We investigated how the activity fraction varies as a func¬ 
tion of rotation period by dividing the sample into four 
logarithmically spaced rotation bins (<1 days, 1-10 days, 
10-100 days, and >100 days). We also used the Adap¬ 
tive Kernel Density Estimation routine akj within the 
quantreg package in R to construct probability density 
functions (PDF) for the a ctive and inactive stars as a 
funct ion of rotation period (jR, Core Teaml[20T3 iKoenkeH 
1201511 . We then combined the active and inactive PDFs 
to compute a PDF for the activity fraction as a function 


We used the Jeffreys prior, 1/a, also known as the logarithmic 
prior, for the standard deviation and a fiat prior for the mean. 


of rotation period. We determined confidence intervals 
(68% and 95%) by taking 5000 bootstrap samples of the 
data. Figure [5] shows the activity fraction as a function 
of rotation period for the logarithmic bins (diamonds) 
and the nonparametric PDF (gray shading). Figure [B] 
indicates that all M dwarfs with periods shorter than 
ten days show detectable Ha emission, with a decrease 
in activity fraction as a function of increasing rotation 
period. While the slowest rotating M dwarfs have lower 
activity fractions, there appears to be a small population 
of active, slowly rotating stars. 



Fig. 5.— The fraction of magnetically active stars for all M dwarfs 
in our sample as a function of rotation period. Diamonds represent 
logarithmically-spaced bins with error bars calculated from bino¬ 
mial statistics. The numbers above/below each symbol represent 
the number of stars in each bin. The gray shaded regions show the 
nonparametric probability density function of the activity fraction 
with confidence intervals (dark gray - 68%; light gray - 95%) de¬ 
termined from 5000 bootstrap samples of the data. Rug plots on 
the top and bottom indicate the active and inactive (respectively) 
stars used to compute the activity fraction. All of the fast rota¬ 
tors are active, whereas a small but non-zero fraction of the slow 
rotators show magnetic activity. 


We further explored the rotation dependence on activ- 
ity fraction by dividing the sample into early (M1-M4) 
and late-type (M5-M8) samples. Figure [B] was made us¬ 
ing the same procedures as in Figure [5] except that we 
used 2000 bootstrap samples to determine the PDF confi¬ 
dence intervals. Figure |6] demonstrates a clear difference 
in the populations with the early-type M dwarfs (left) 
showing a strong decrease in activity fraction with in¬ 
creasing rotation period, with the slowest rotators being 
completely inactive. The activity fraction of the late- 
type population (right) stays very high for stars with 
rotation periods shorter than 100 days, after which it de¬ 
creases. A small number of the slowest rotating, late-type 
M dwarfs still show magnetic activity. While differences 
in the activity fractions between early and late-type M 
dwarfs have been seen for decades fe.g.. lGizis et al.ll2000l : 
iWest et al.ll200^ and can be exp lained by differences in 
active lifetimes (jWest et'ani2008h . Figure |6] suggests that 
the relationship between activity and rotation may be 
different in late-type versus early-type M dwarfs. Specif¬ 
ically, it appears that late-type M dwarfs may stay active 
at slower rotation speeds than their early-type counter¬ 
parts. 
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M1 - M4 M5 - M8 



0.1 1.0 10.0 100.0 0.1 1.0 10.0 100.0 
Period (days) Period (days) 


Fig. 6.— The fraction of active early-type (left) and late-type (right) M dwarfs in our sample as a function of rotation period. Diamonds 
represent logarithmically-spaced bins with error bars calculated from binomial statistics. The numbers above/below each symbol represent 
the number of stars in each bin. The gray shaded regions show the nonparametric probability density function of the activity fraction 
with confidence intervals (dark gray - 68%; light gray - 95%) determined from 2000 bootstrap samples of the data. Rug plots on the top 
and bottom indicate the active and inactive (respectively) stars used to compute the activity fraction. At similar rotation periods, a much 
larger fraction of the late-type M dwarfs are active, evidencing that the activity-rotation relation in M dwarfs may be mass dependent. 


We also investigated how the strength of magnetic ac¬ 
tivity (quantified by Ana/Aboi) varies as a function of 
rotation period. Figure[7]shows log(LHa/Aboi) (filled cir¬ 
cles) as a function of rotation period for early-type (left; 
M1-M4) and late-type (right; M5-M8) M dwarfs. The in¬ 
active stars with measured rotation periods are plotted as 
open circles. The early-type M dwarfs show a decrease in 
activity strength with increased rotation period. In con¬ 
trast, the late-type M dwarfs are consistent with having 
the same level of activity, except perhaps at the longest 
rotation periods. To quantify this, we performed a lin¬ 
ear least squares fit to the active stars in both panels of 
Figure [T] In the early-type dwarfs, log(THa/Tboi) shows 
a statistically significant decrease as function of rotation 
period, and has the form, 

log(LHa/Tboi) =-3.44 ±0.02 

-0.19 ± 0.04 xlog(Prot), (3) 

with Prot as the rotation period in days. For the late- 
type M dwarfs, log(LHa/Tboi) shows no significant trend 
with rotation period (slope = —0.016 ±0.050). The slope 
of the active late-type dwarfs is not consistent with the 
decreasing trend seen in the active early-type M dwarfs. 
In both cases, the scatter about the best-fit line is dom¬ 
inated by the intrinsic scatter rather than measurement 
uncertainty. Binning the data in logarithmic bins of rota¬ 
tion period confirms these bulk trends (Figure [8]). Both 


Figures [7] and [8] demonstrate that there is a connection 
between rotation and activity in late-type M dwarfs; all 
stars rotating with periods shorter than ^90 days are 
active. There is also a hint that there may be a corre¬ 
lation between rotation and activity in the late-type M 
dwarfs with rotation periods longer than 90 days, but 
there are not enough stars to make any strong conclu¬ 
sions (e.g., the slowest rotation bin of Figure |8] contains 
only one star). Additional late-type, slowly rotating, ac¬ 
tive M dwarfs are required to confirm if magnetic activity 
is weaker in slowly rotating, fully convective M dwarfs. 

Figure [7] also suggests that while fast rotation is cor¬ 
related with magnetic activity in all M dwarfs, there are 
clear differences between the early and late-type popu¬ 
lations. Similar to what was seen in Figure [HI the late- 
type M dwarfs maintain a high (and comparable) level of 
magnetic activity until their rotation periods exceed ^90 
days, beyond which the level of activity appears to de¬ 
crease. In the early-type M dwarfs, there is a much more 
clear decrease in activity at all rotation periods and an 
absence of active stars with periods > 30 days. The inac¬ 
tive stars in both panels corroborate that rotation plays 
an important role in the generation of magnetic activ¬ 
ity in M dwarfs and hint at the possibility of a rotation 
threshold faster than which magnetic activity is present. 
All M dwarfs rotating faster than 26 days and 86 days 
in the early-type and late-type populations respectively 
are magnetically active. 
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M1-M4 M5-M8 



Fig. 7.— log(LHa/tjbol) S'S a function of rotation period for early-type (left) and late-type (right) M dwarfs (filled circles). Inactive 
stars with measured rotation periods are included as open circles. In the early-type M dwarfs there is a clear decease in the strength of 
magnetic activity with increasing rotation period. The late-type M dwarfs are much more scattered, but appear to maintain a similar level 
of activity except for the very slowest rotators. The inactive stars suggest that there may exist rotation periods (26 days; early-type and 
86 days; late-type), faster than which all stars are magnetically active. 


M1-M4 M5-M8 



Rotation Period (days) Rotation Period (days) 

Fig. 8.— log(THa/Tboi) a function of rotation period for early-type (left) and late-type (right) M dwarfs. Each data point corresponds 
to the mean value (of activity and rotation period) for the stars in the bin. The vertical error bars represent the spread of the data (thin 
bars) and the uncertainty in the mean values (wide bars). Only three bins are shown for the M1-M4 dwarfs, as no active stars of these 
types were found with periods longer than 100 days. There is only one star in the slowest rotating late-type bin and the error bars on 
this data point represent the uncertainty in the I/HaZ-^bol measurement. In the early-type M dwarfs there is evidence of a decease in the 
strength of magnetic activity with increasing rotation periods. This is in contrast to the late-type M dwarfs, which maintain a similar level 
of activity except for the very slowest rotators. 
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Fig. 9.— The mean V velocities (top), and W velocity dispersions 
(bottom) as a function of rotation period. Error bars represent 
the IfT uncertainties in each bin. Both panels suggest that slowly 
rotating M dwarfs are drawn from an older kinematic population. 

Figures [7] and [8] also demonstrate an alternative 
method for studying the range of magnetic activity 
strengths in M dwarfs. Our results reduce the contribu¬ 
tion of different rotation periods to the spread of activity 
strength and allowed us to examine the intrinsic spread 
of the population (at a given rotation rate), particularly 
for early-type M dwarfs. While the mean values have rel¬ 
atively small uncertainties (wide error bars) in most bins, 
the spread (thin error bars) at each rotation bin remains 
large (typically about 0.3-0.5 dex). We see that mix¬ 
ing stars of various rotation rates (and perhaps different 
ages) does increase the spread of the population (particu¬ 
larly for early-type M dwarfs). Additionally, some of the 
spread seen in individual bins in Figure [8] indicate that 
the M dwarf population has a significant intrinsic variety 
of activit y strengths. As has been seen in previous stud¬ 
ies (e .g.. iBerger et al. I [20081 : iLee et all 120101 : iBell et aP 
I2012D . much of the observed spread may be due to the 
Ha variability of individual M dwarfs. 

4.2. Kinematics, Rotation and Activity 

From the results of the kinematic analysis described in 
Section 3.2, we examined the bulk motions of the sample 
as a function of rotation period. Figure[9]shows the mean 
V velocity (top) and the W velocity dispersion (bottom) 
of the entire sample (relative to the Local Standard of 
Rest) as a function of rotation period. The V compo¬ 
nent of the velocity (in the direction of Galactic rotation) 
shows clear signs of increasing asymmetric drift with in¬ 
creasing rotation period. Namely, the mean V velocity 
decreases for stars with slow rotation rates. This de¬ 
crease is a strong indicator that the population of slow 
rotators is older; dynamical heating over time causes el¬ 
liptical orbits, which lead to a bulk decrease in V motion. 


The mean U and W motions are insensitive to asymmet¬ 
ric drift and provide little information about the ages of 
small samples of stars. All of the population U, V, and 
W values (means and dispersions) are included in Table 
[2] for completeness. 

The velocity dispersion of the W component (bottom 
panel of Figure IHl in the vertical direction, out of the 
plane of the Galaxy) increases with increasing stellar ro¬ 
tation period. Increases in the vertical velocity disper¬ 
sion are typically associated with multiple dynamical en¬ 
counters and older ages. As with the asymmetric drift 
apparent in the mean velocities, the vertical velocity dis¬ 
persion measurements suggest that the slower rotators 
are dynamically older. The dispersions of the U and V 
components can be affected by other processes and are 
less directly linked to dynamical heating than the W ve¬ 
locity dispersions, but are also included in Table [2j 

We also completed the same kinematic analysis for sep¬ 
arated late-type and early-type populations and found no 
significant differences in their kinematics. 

5. CONCLUSIONS AND DISCUSSION 

The spectroscopic follow-up of 238 M dwarfs, combined 
with the rotation rates from the MEarth transit survey 
has provided an unprecedented sample for examining the 
rotation, magnetic activity and kinematics of mid-to-late 
M dwarfs. We conducted an analysis of the magnetic ac¬ 
tivity properties and investigated the bulk 3D kinemat¬ 
ics of the population as a function of rotation. In the 
spectroscopic sample, 164 stars have measured rotation 
periods, 160 are magnetically active, and 127 of the M 
dwarfs are both active and have measured rotation peri¬ 
ods. We provide all of the sample and derived quantities 
in the electronic version of this manuscript. 

The results of our analysis are summarized as follows: 

1. For early-type M dwarfs, the activity fractions de¬ 
crease with increasing rotation period. For the late- 
type M dwarfs, the higher activity fractions extend 
to slower rotation rates before rapidly declining. 

2. All of the rapidly rotating stars are magnetically 
active, with few (but some) of the slowest rotators 
showing signs of activity. 

3. For early-type (M1-M4) M dwarfs, all stars rotat¬ 
ing faster than 26 days are magnetically active, as 
are all late-type M dwarfs (M5-M8) with rotation 
periods shorter than 86 days. It is worth noting 
that for a 0.2 Mq star, 26 and 86 days correspond 
to V sin i values of 0.4 and 0.1 km s“^ respectively, 
both of which are below the detection limits of most 
spectrographs. 

4. At the longest rotation periods, we see no active 
early-type M dwarfs, but there are a few slowly 
rotating, active, late-type M dwarfs. 

5. The strength of magnetic activity appears to de¬ 
cline with increased rotation period in early-type 
M dwarfs. The fully-convective, late-type M dwarfs 
remain at a similar activity level as rotation slows, 
perhaps until the slowest rotation periods (> 90 
days). 
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6. From our kinematic analysis, it appears that M 
dwarfs with longer rotation periods are drawn from 
a kinematically older population. Specifically, we 
see that slowly rotating stars exhibit signs of asym¬ 
metric drift and have larger velocity dispersions 
about the mid-plane, both signs of dynamically 
heated, older populations. 

7. There are a few slowly rotating M dwarfs that show 
signs of magnetic activity, and which may provide 
important laboratories for additional investigations 
about the role rotation plays in the magnetic field 
generation in low-mass stars. 

In general, we find a strong relationship between rota¬ 
tion and activity for early-type M dwarfs. The relation¬ 
ship is much less pronounced, but appears to exist for the 
late-type M dwarfs that are fully convective. For all M 
dwarfs, rotation appears to affect the presence of activ¬ 
ity; fast rotators are more likely to show strong activity, 
while slow rotators are less likely to be active. However, 
in the late-type M dwarfs, magnetic activity appears to 
be present in stars with slower rotation than their early- 
type counterparts. Larger samples that include measure¬ 
ments of magnetic activity and rotation are required to 
further investigate this discrepancy (e.g. Newton et al. 
in prep). 

Our kinematic analysis suggests that rotation is linked 
to age in all M dwarfs. We conclude that younger M 
dwarfs (both early and late-type) are rotating faster and 
more likely to be active. Future analysis of wide binaries 
(particularly those with a white dwarf companion from 
which we can determine age) and/or stellar clusters will 
provide better information about the detailed role that 
stellar age plays in the observed rotation and magnetic 
activity of low-mass stars. 

While our kinematic analysis shows clear signs of asym¬ 
metric drift and dynamical heating, the magnitudes of 
the effects are far less than what is seen in large sam¬ 
ples of M dwarfs that probe farthe r into the Milky Way 
disk (e.g., iBochanski et al.lF2007all . This suggests that 
the slow rotators in our sample are older than the fast 
rotators, but that they are still younger than typical old 
stellar populations in the Galactic (thin) disk. This is 
not surprising given that the MEarth targets are selected 
to be nearby, where it is much more likely to sample a 
younger population (jWest et al.ll2008[) . 

Some of the active stars in our spectroscopic sample 
do not have measured rotation periods (21%) and were 
therefore not included in some of our analyses. 74 stars 
in the sample do not have rotation periods (55% of which 
are inactive). Because our ability to measure a rotation 
rate of an M dwarf depends on the presence of starspots, 
it is not surprising that some stars, namely the inac¬ 
tive ones, may not show rotational modulation; indeed, 
the majority of the stars with rotation periods are active 
(77%). However, the lack of measured rotation in the ac¬ 
tive stars could be due to a number of effects, including 
but not limited to, homogeneous spot geometry, short¬ 
lived spots, insensitivity to low-amplitude variations (few 
small spots), or an over abundance of spots. The latter 
effect would cause a decrease in the rotational modula¬ 
tion for the most active (heavily spotted) stars. We ran 
the non-rotating stars through our kinematic analysis (53 


have 3D kinematics) and found that as a population they 
are consistent with M dwarfs in the 10-100 day rotation 
bins. We conclude that kinematically, the stars without 
rotation measurements could be on the older end of the 
spectrum for our sample. The active fractions of early 
and late-type M dwarfs without rotations measurements 
are consistent with the entire sample in the 10-100 day 
rotation bins. Therefore, the exclusion of the active stars 
without rotation measurements does not likely effect our 
results, but may be explored in a future study when we 
can ascertain the reason for the lack of rotational modu¬ 
lation. 

The fact that there are a small number of active, slowly 
rotating, late-type M dwarfs in our sample suggests that 
the activity-rotation relation may be more complicated 
than predicted by a simple spin-down model. One pos¬ 
sibility could be that while persistent magnetic activity 
in late-type M dwarfs is rare (or non-existent) for slow 
rotators, the presence of magnetic cycles may produce 
observed activity in a small number of stars at any given 
time. Long-term observations are required to further test 
this hypothesis. The small number of active, slowly ro¬ 
tating, late-type M dwarfs could also be due to the young 
bias of MEarth, which might exclude such stars. An ad¬ 
ditional factor that could influence our analysis is the po¬ 
tential presence of close binary companions, which have 
been shown to affe ct the magnet i c act ivity (and rota¬ 
tion) of stars le.g.. [Morgan et all 1201211 . While multi¬ 
plicity likely plays a minor role due to the sm all binary 
fraction of M dwarfs ([Fischer fc Marc^ll992ll . high res¬ 
olution imaging and/or spectroscopy would be required 
to further investigate its role. 

Our analysis divided M dwarfs into early (M1-M4) 
and late-types (M5-M8) and did not explore the ram- 
ifcations of altering the boundary between the two pop¬ 
ulations. Ideally, we would conduct the analysis as a 
function of spectral type but we did not have enough 
stars in each spectral type to produce a statistically sig¬ 
nificant result. As stated above, our analysis is there¬ 
fore mostly a comparison between M3-M4 and M5-M6 
dwarfs. Future and ongoing studies will provide addi¬ 
tional time-series photometric data with which we can 
further study the rotation and activity behavior of M 
dwarfs (over a larger range of spectral types ), both sur¬ 
veys such as MEart h (both North and So uth: llrwin et ahl 
[ml, APACHE ISozzetti et af.l [2?fn . SPECULOOS 
([Gillon et al.1 I2013I1 . and ExTrA[3 that focus specifi¬ 
cally on M dwarfs and more general photo metric sur- 
yeys like Pan-STA RRS ([Kaiser et al.l l2002tl and LSST 
([Ivezic et al.ll200l . 
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TABLE 1 

Kinematics, Activity, and Rotation Periods for MEarth M Dwarfs'^ 


Name 

RA 

Dec 

Sp. 

Type 

Dist^^ 

(pc) 

Dist^^ 

Flag 

RV 

(km/s) 

PMra‘> 

(mas/yr) 

PMDec'" 

(mas/yr) 

(km/s) 

(km/s) 

(km/s) 

Per. 

(days) 

Ha EW'i 

(A) 

Act. 

Flag^ 

Rot. 

Flagf 

LSPMJOOOl+0659 

00;01:15.809 

+06:59:35.65 

M6 

16.7 

D 

-0.7 

-447 

-81 

+11.1 ± 1.5 

±0.1 ±8.4 

±13.8 + 4.6 

20.4 

±6.9 ±1.2 

1 

1 

LSPMJ0015+4344 

00:15:18.827 

+43:44:34.72 

M5 

25.7 

L 

+5.9 

+232 

+38 

+10.1 + 4.1 

+4.6 ±8.3 

+7.1 ±2.8 

1.4 

±6.5 ± 0.8 

1 

1 

LSPMJ0015+1333 

00:15:49.244 

+13:33:22.25 

M3 

9.6 

D 

+65.1 

+621 

+333 

+0.2 ±2.2 

±32.1 + 8.3 

-25.4 + 4.3 



0 

0 

LSPMJ0016+1951E 

00:16:16.141 

+19:51:50.56 

M4 

22.2 

L 

+9.0 

+709 

-748 

+10.3 ±2.8 

+4.1 ±8.3 

+7.9 ±4.1 



0 

0 

LSPMJ0016+2003 

00:16:56.803 

+20:03:55.17 

M4 

22.2 

L 

+12.0 

+228 

+24 

±9.2 ±2.8 

+6.9 ±8.3 

+5.2 ±4.0 

17.3 

±3.5 ± 0.4 

1 

1 

LSPMJ0018+2748 

00:18:53.590 

+27:48:49.81 

M4 

19.2 

L 

+12.7 

+387 

-101 

±8.5 ±3.3 

+8.3 ±8.3 

+4.5 ±3.6 

6.0 

±1.6 ± 0.3 

1 

1 

LSPMJ0024+2626 

00:24:03.799 

+26:26:29.76 

M4 

26.1 

D 

+16.0 

+162 

-55 

±8.3 ±3.5 

+8.5 ±8.3 

+4.2 ±3.7 

29.9 

±1.3 ± 0.4 

1 

1 

LSPMJ0024+3002 

00:24:34.876 

+30:02:29.59 

M5 

15.5 

D 

+17.6 

+580 

+28 

±7.1 ±3.6 

±10.8 + 8.3 

+2.9 ±3.4 

1.1 

±6.6 ± 0.9 

1 

1 

LSPMJ0028+5022 

00:28:53.972 

+50:22:33.17 

M4 

16.9 

D 

+24.2 

+423 

+124 

±2.8 ±4.6 

±17.8 + 8.3 

+3.8 ±1.8 

1.1 

±5.9 ± 0.6 

1 

1 

LSPMJ0033+1448 

00:33:22.350 

+14:48:06.45 

M5 

27.3 

D 

+21.5 

+278 

-71 

±6.7 ±3.1 

±11.6 + 8.2 

-1.2 ±4.2 

0.4 

±3.5 ± 0.6 

1 

.2 

LSPMJ0035+5241N 

00:35:53.680 

+52:41:36.59 

M4 


L 

+12.7 

+787 

-186 






0 

§ 0 

LSPMJ0038+6150 

00:38:27.677 

+61:50:06.34 

M4 

46.5 

D 

-58.3 

+349 

-43 

±37.5 + 5.1 

-40.1 + 8.2 

+7.9 ±1.2 

62.4 


0 

fD o 

CO ^ 

LSPMJ0039+1454S 

00:39:33.544 

+14:54:18.96 

M4 


L 

+9.1 

+321 

+39 




34.0 


0 

^ 1 

LSPMJ0103+6221 

01:03:19.824 

+62:21:55.74 

M6 

10.5 

L 

+6.7 

+739 

+86 

±9.3 ±5.4 

+5.9 ±7.9 

+7.0 ±0.5 

1.0 

±9.2 ± 1.3 

1 

S. 1 

LSPMJ0130+0236 

01:30:43.122 

+02:36:37.01 

M6 

18.3 

D 

-15.6 

-119 

-194 

±18.3 + 4.1 

-1.5 ±7.6 

±24.3 + 4.6 


±0.8 ± 0.5 

1 

0 

LSPMJ0153+0147 

01:53:30.754 

+01:47:55.89 

M6 

20.1 

D 

+30.4 

+427 

+45 

±1.4 ±4.7 

+9.4 ±7.2 

-8.2 ±4.6 

0.2 

±8.6 ± 1.2 

1 

■ 1 

LSPMJ0153+4427 

01:53:49.551 

+44:27:28.49 

M5 

20.1 

D 

+14.3 

+245 

-92 

±6.4 ±6.5 

+8.6 ±7.1 

+5.4 ±0.7 

0.2 

+17.8 ± 2.3 

1 

1 

LSPMJ0158+4049 

01:58:45.211 

+40:49:44.51 

M5 

17.1 

D 

+19.1 

+399 

-69 

±3.5 ±6.5 

±11.1 ± 7.0 

+3.7 ±0.9 

0.5 

+7.3 ± 1.0 

1 

1 

LSPMJ0202+1334 

02:02:44.355 

+13:34:33.21 

M5 

20.7 

D 

+15.6 

+461 

-111 

±6.7 ±5.7 

+7.0 ±7.0 

+3.0 ±3.5 

4.0 

+7.6 ± 0.9 

1 

1 

LSPMJ0212+0000 

02:12:54.624 

+00:00:16.71 

M4 

9.8 

D 

+33.0 

+560 

+34 

±0.5 ±5.1 

+8.0 ±6.7 

-8.1 ±4.7 

4.7 

+2.2 ± 0.4 

1 

1 


- This table is published in its entirety i 


the electronic edition of the Astrophysical Journal. A portion is shown here for guidance concerning its form and content. The machine-readable table also contains the semiamplitudes 

et al.l 


of the rotation data as well as other extracted parameters from the spectra used in this paper, including a calculation of and atomic/molecular indices for Na, TiO, CaH, CaOH features (as defined in 

^ A Distance Flag of ”D” indicates the distance is a trigonometric parallax measurement derived from MEarth imaging, as published in[Dittman^_et_a^ . A flag of ”L” indicates the distance comes from 

may be either photometric, spectroscopic, or trigonometric. 

Proper motions are quoted as projected on the plane of the sky, with (PMj^^, PMj^g^.) = (^q.cos 5, p,^). 

^ UVW velocities are quoted in a right-handed coordinate system, with U pointed towards the Galactic center. 

^ This table uses a convention in which emission lines correspond to equivalent widths > 0. 

^ Active M dwarfs have an Active Flag value of 1; inactive stars have a flag of 0 (see Section 

Stars with measured rotation periods have a Rotation Flag value of 1 or 2, with a value of 1 being more robust. Stars whose rotation periods did not cross our detection threshold are flagged as 0 (see Section[2^^. 


I i2QQai and 
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TABLE 2 

Population Kinematics - Entire Sample 


Mean Rotation 
Period (days) 

f^U 

(km s“^) 

(km s“^) 

fiW 

(km s“^) 

(km s“^) 

CTV 

(km s“^) 

(km s“^) 

0.5 

3.3 

49.0 

119.8 

9-8JI 

o qO.7 

K fjl.O 
92!:t 

144^ 

8-8M 

14.2Jg 

6.3?:! 

4 c^O.6 

84!:e 

<19. 

'-^1.7 

7 oO.S 

V-? 

9.8il 

14.41 






